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We consider the possible role of rescattering effects in the determination of the gluon density 
for the LHC from DIS data. We discuss a method that uses results of s-channel calculations to 
estimate these effects, and comment on potential applications to diffractive and multi-parton 
interactions. 



The Large Hadron Collider will operate with very high gluon luminosities. Production 
processes initiated by gluons contribute a great many events to a number of cross sections of 
primary interest for the LHC physics program. Reliable predictions for these cr oss secti ons 
depend on the determination of the gluon density in the proton and its accuracy. HEEEE] ^ g 
parton luminosities rise steeply for decreasing momentum fraction x, a large number of events 
i—i. sample the gluon density at x <C 1. 

'"q! ^ nas long been known that the theoretical accuracy of gluon-density determinations for 

x <C 1 is affected by higher-loop ln(l/x) corrections to QCD evolution equations. See e.g. Ref.^ 
for early numerical investiga tions. The study of these corrections motivates current searches for 
evolution schemes (see Refs. 1 1 1 \ and references therein) that i ncorpo rate the resummation of 
ln(l/x) contributions at the next-to-leading-logarithmic accuracy. An improved theoretical 

r—i ■ control on the xCl region is expected from the inclusion of these terms. 

Because the DIS data used at present to extract the gluon density for x < 10~ 2 do not have 

£f~) | very high Q 2 , it is natural to ask whether non-negligible effects on the theoretical accuracy may 

also come from corrections that are suppressed by powers of 1/Q 2 but are potentially enhanced 
as x — > 0. These could affect the determination of the gluon density f g primarily through a 
contribution 5 to the Q 2 -derivative of the F% structure function, 

■ dp 1 a 2 
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Here P qg is the perturbative gluon-to-quark evolution kernel, and the correction 5 arises from 
multi-parton correlation terms in the operator product expansion, 

F 2 = C®f + ^C^®f^ + ... . (2) 

The enhanced x — > behavior in Eq. ([T]) can be produced from graphs with multiple gluon 
scatterings, and is consistent wi th observa tions of approximate geometric scaling in low-x data^l 
and with models for saturation J 13 | 14 115| 

Standard methods to take account of multiple scatterings are s-channel methods (see e.g. 
the lectures in Ref.^1), essentially orthogonal to those of the parton picture. The basic degrees 
of freedom in the s-channel picture are described by correlators of eikonal Wilson lines — at the 
simplest level, two-point correlators, interpretable as color dipoles. To identify the correction 
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from rescattering graphs to the parton result in Eq. ([T]), a sufficiently precise "dictionary" is 
needed to connect the two pictures. Ref.^l presents an approach to analyze this connection. 

The method is based on constructing explicitly an s-channel representation for the renormal- 
ized parton distribution function in terms of Wilson-line matrix elements, convoluted with light- 
cone wave functions. In this representation the quark distribution f q is given by the coordinate- 
space convolution 

xf q (x,n) = J dz J db u(n,z) E(z,b) , (3) 
where 5 is the hadronic matrix element of eikonal-line operators, 

E(z,b) = J[dP'] (P'\±-Tr{l-V\b + z/2)V(b-z/2)}\P) , 

r r+oo ^ (4) 

V(z) =Vexpl-igJ dz-A+{Q,z-,z)t a \ , 

z is the transverse separation between the eikonal lines, b is the impact parameter, and the 
function it(/i, z) is evaluated explicitly in Ref. ^1 a t one loop using the MS scheme for the 
renormalization of the ultraviolet divergences z — > 0. 

The representation ([3]), once ev aluated in a well-prescribed renormalization scheme, is the key 
ingredient that allows one to relateMZl 

results of s-channel calculations for structure functions 
to the OPE factorization (J2j) , and, in particular, to identify the power-suppressed corrections 
arising from the multiple gluon scatterings, treated in the high-energy approximation of Eq. (JH). 
These corrections are found to depend on moments of E, schematically in the form 

dF 2 ( dF 2 \ ~ A 2 (n) 



dlnQ 2 \dlnQ 2 J Lp ^ (Q 



where the first term in the right hand side is the leading-power parton result, and the A 2 in the 
subleading terms are given by the analytically continued moments 

A 2 (-) = ^ / ^ (^r 1 JdbS(z,b) . (6) 

The coefficients R n are evaluated to order a s , as functions of x and InQ 2 , from the lightcone 
wave functions, while the moments A 2 (n) are dimensionful nonperturbative parameters, to be 
determined from comparison with experimental data. 

In practice, the usefulness of the result in Eqs. ©,© comes from the fact that the hadronic 
matrix element 5 can be related by a short-distance expansion for z — > to a well-prescribed 
integral of the gluon distribution functional. Then the moments A 2 can be parameterized in 
terms of the factorization/renormalization scales at which the gluon distribution and the strong 
coupling are evaluated. These scales are to be taken of the order of the inverse mean transverse 
distance l/\z\, and can be tuned to the data. The result of doing this for F 2 data^^at both low 
and high Q 2 is shown in Fig. [T] in the left hand side plot.Q^ The corresponding power correction 
in Eq. §5§ is plotted on the right hand side of Fig. [TJ Here we show the correction normalized 
to the full answer and multiplied by —1, using CTEQ parton distributions.^^ 

Fig. Q] indicates that with physically natural choices of the parameters in the nonperturbative 
matrix elements (jlj) one can achieve a sensible description of data for x < 1CP 2 in a wide range 
of Q 2 and still have moderate power corrections to dF 2 /dlnQ 2 . Corrections turn out to be 
negative and below 20 % for x > 10 -4 and Q 2 > 1 GeV 2 . However, Fig. Q] also indicates that for 
small x the corrections fall off slowly in the region of medium Q 2 , Q 2 ~ 1 — 10 GeV 2 , behaving 
effectively like 1/Q X with A close to 1. ^ For instance, one has A ~ 1.2 for the curve x = 10~ 3 
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Figure 1: (left) The result of fitting the A 2 parameters to the data [18]; (right) power corrections to dF^/dlaQ 2 

versus Q 2 at different values of x. 



The rationale for this expansion lies with the dynamical cut-off on large distances z imposed 

I H J I -J M l -j A\ 

by unitarity requirements ("black disc" limit) on the correlator 5j ' ' ' But the size of this 



in the right-hand side plot of Fig. [TJ As a consequence of the slowly decreasing behavior, the 
power corrections stay on the order of 10% up to Q 2 of a few GeV 2 for x < 10~ 3 . This slow 
fall-off differs from parameterizations of higher twist commonly used in global analyses (see e.g. 
Ref.^), and may be relevant for phenomenology as it affects the medium Q 2 region of the data 
used to extract f g at low x. 

Corrections larger than for F% are found for the longitudinal component Fl This provides 
additional motivation for the forthcoming Fl measurements^^, as well as fits^H investigating 
power- like terms in Fl (at both high and low x). We observe also that Fig. Q]is obtained using 
NLO parton distributions, and the decrease in the low-x gluon at NNLO could be consistent 
with the possibility that NNLO parton distributions correspond to smaller power corrections. 
However, the detailed interpretation of this behavior will be subtle, as distinctly dif feren t dy- 
namics drive the power-like and NNLO effects, unlike the high-x case in the analyses 15^1 . 

It is worth emphasizing that the above results depend on the validity of the high-energy 

approximation and s-channel representation ([3|), and the perturbation expansion for u(fi,z). 

T on large d 
5 J14ll5Tl6l 

cut-off at collider energies is difficult to determine. The highest se nsitivit y to it may come 

91 24 9% 

from measurements of the diffractive part of the DIS cross section. 1 111 In this case the 
s-channel representation is bilinear in S. The comparison^^Sl of diffractive data with theoretical 
predictions based on diffractive parton distributions indicates that the dynamical cut-off lies at 
substantially higher momenta for color-octet eikonal-line matrix elements than for color-triplet. 
That is, gluons' shad ow is s tronger. See e.g. Ref. for a recent discussion. In diffractive 
DIS this can be linked 1 ^™ to the distinctive pat tern of th e observed scaling violation^23! an( j 
detailed features of the associated jet distributions. ^SEHEQ] More generally, it suggests that the 
expansion used is better justified for processes directly coupled to the gluon distribution than 
for F2, see e.g. applications to Fl (or its diffractive component^^) and jet final states. A critical 
discussion, including the quark case, is given in Ref.^^. 

Note that the question of how to perform QCD calculations that incorporate multiple scat- 
terings along with perturbative evolution becomes especially compelling in the case of Monte 
Carlo event generators^]. The application considered above deals with corrections to Q 2 evo- 
lution, i.e., a picture based on strongly ordered k± , s. But the method relies on high-energy 
approximations that may be better suited for extension to evolution with ordering in energies, 
or angles. It could thu s mor e likely be adapted to the modeling of multiparton processes in 
Monte Carlo generators 1 1 1 based on high-energy evolution equations. 
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